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ABSTRACT 

There are indications that the neutrino oscillation data from reactor experiments and the LSND and MiniBooNE experiments show 
a preference for two sterile neutrino species, both with masses in the eV region. We show that this result has a significant impact 
on some important cosmological parameters. Specifically, we use a combination of CMB, LSS and SN1A data and show that the 
existence of two light, sterile neutrinos would rule out the cosmological constant as dark energy at 95% confidence level, and lower 
the expansion age of the universe to 12.58 ± 0.26 Gyr. 
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1. Introduction 

It is no news that the properties of neutrinos impact the history 
of the universe. For example, the first limits on the number of 
neutrino flavours cam e from Big Bang Nucleosynthesis (BBN) 
dSteigman et al.llT977l) . And as measurements of the statistical 
distribution of matter in the universe and the anisotropics in the 
cosmic microwave background (CMB) have improved, it has be- 
come possible to put increasi ngly stringent upper bounds on the 
sum o f the neutrino masses (Ko matsu et alJ[20TTt iThomas et alJ 
l2010h . The strongest bounds result, of course, when one starts 
from the simplest cosmological model with a handful of param- 
eters fitted to a selection of the most important data sets and then 
includes the sum of the neutrino masses as an additional degree 
of freedom. 

However, there are good reasons to take the cosmological 
mass limits with a grain of salt. For one thing, the cause of 
the apparent accelerated expansion of the universe is unknown. 
Although th e cosmologica l cons tant is consistent with all ex- 
isting data dKomatsu et a n i2oTTh . we cannot exclude alterna- 
tive explanat ions like scalar fi elds dMota & Shaw 2007) , modi - 
fied gravity (Tsujikawa 2010) or void models (lMattssonll2010l) . 
Furthermore, neutrino experiments suggest that the neutrino sec- 
tor of the Standard Model of particle physics may be more com- 
plicated and interesting than the simplest picture with three mas- 
sive flavour s tates consisting of a superposition of three mass 
eigenstates (lAguilar et al.l 120011: lAguilar-Arevalo etal 1 I201C1 
Muel ler et al.l 1201 H) . In this paper our point of departure is 
the fact that the simplest realizations of massive neutrinos do 
not seem to explain the results of experiments with neutrinos 
from nuclear reactors and the LSND and MiniBooNE exper- 
iments. Recently it was suggested that two sterile neutrinos, 
neutrinos that only interact gravitationally with matter, give a 
ood description of the data if their masses are in the eV-range 
Koppet alhoilh . 

The ability of present and future cosmological data 
sets to constrain light sterile neutrinos was i nvesti gated 
in (lGonzalez-Garcia etal.ll2010h . ICalabrese et al.l (l201ll) . and 
Giusarma et al.l (1201 ll) . We choose to approach the problem from 
the opposite direction: given the uncertain nature of the dark en- 
ergy, and the somewhat model-dependent interpretation of cos- 



mological data, we would argue, like we have done in the past 
dKristiansen & Elgar0vll2OO8l: iKristiansen etaT]l2010h . that neu- 
trino experiments have greater authority than cosmology. If the 
latter find that light, sterile neutrinos are required, cosmolo- 
gist have to find room for them in their models. Knowing full 
well that the question is far from settled, we nevertheless find it 
worthwhile to consider what changes in the cosmological con- 
cordance model a scenario with 3 active and 2 light, sterile neu- 
trinos lead to. We will therefore i nvestigate the scenario where 
the proposal in ( Kopp et alJl20lTl) is assumed to be correct and 
factor their result into an analysis of current cosmological data. 

Our paper is organized as follows. In section II and III we 
summarize the theoretical background and describe our method, 
but briefly since it follows procedures that are standard in the 
literature. Section IV is the most important section where we 
present our results and the inferences we draw from them. We 
summarize and conclude in section V. 

2. Sterile neutrinos and cosmology 

The number of neutrino species with masses below the 
GeV scale and that couple to the Z° boson, i.e., interact 
weakly, was determined to be 2.984 + 008 from LEP data 
dNakamura & Particle Data G roup 2010). If there are more neu- 
trino types than the three we already know about, they must be 
very heavy, or couple to gravity only, or both. Neutrinos that do 
not participate in the weak interaction are k nown as ste rile. They 
appear in the so-called seesaw mechanism (Zuber 2004) for gen- 
erating small neutrino masses, and are there typically very heavy, 
much heavier than the electroweak scale, in order to explain the 
smallness of the masses of the ordinary, active neutrinos. 

However, as long as it only interacts gravitationally there 
are no a priori constraints on the mass of a putative sterile neu- 
trino. Sterile neutrinos with keV masses have b een of great in- 
terest as dark matter candidates in cosmology dKusenkol F2009) 
. And recently it has been suggested that one or two sterilie 
neutrinos with masses of a few eV lie beh ind some puzzlin g 
features in neutrino oscillation experiments (Kopp et ~al]l20lll) . 
Ever since the LSND experiment fou nd indications of v M -v e 
transitions (Athanassop oulos et al.lll998l) . there have been spec- 
ulations about the existence of a light, sterile neutrino. The 
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MiniBooNE experiment ( Agu ilar-Arevalo et al .2010) provided 
support for the LSND result, but found no eviden ce for oscilla- 
tions in the v^-Vg channel. A recent re-evaluation( Mueller et al. 
1201 ll) of the expected antineutrino flux from nuclear reactors 
hint at neutrino oscillations over length scales of tens to hun- 
dreds of meters. All of these results have been shown to be ac- 
commodated within a model with two eV-mass sterile neutri- 
nos that are quite strongly mixed with electron-type neutrinos. 
A single sterile species is compatible with all the results except 
the negative MiniBooNE result for the v„-V e. channel. We will 
take the best- fit models of (Kopp et al. 201 1) as our point of de- 
parture, investigating both the models with one and two sterile 
neutrino species. 

An important question when we turn to the cosmological im- 
plications of these two scenarios is whether these light, sterile 
neutrinos were thermalized in the early u niverse. We will as- 
sume that they were, since several studies (lHamann et al.ll20Tot 
iMelchiorri et al.ll2009t iKainulainenll 1 990h suggest that this was 
the case for the masses and mixing parameters we consider. This 
means that standard relation between the sum of the neutrino 
masses and their contribution to the cosmic mass density param- 
eter applies. 

Adding two light sterile neutrinos may cause some prob- 
lems with BBN since the increased relativistic energy density 
results in a larger neutron-to-proton ratio, and leads to an in- 
creased He/H mass fracti on. For example, a recent analysis 
dMangano & Serpicol201 ll) found that BBN constrains the num- 
ber of relativistic degrees of freedom to be A^ff < 4.2 at 95 % 
confidence. So the 3+1 model is just within the bounds, and the 
3+2 model is just outside it. If the neutrino oscillation data ends 
up pointing unequivocally to the existence of two light, sterile 
neutrinos, this would mean that the standard BBN scenario has 
to be modified. This is beyond the scope of this paper. 

In contrast to the relatively tight upper bound on A^ff from 
BBN, several recent papers have suggested that additional rel- 
ativistic species are allowed, and in fact preferred, by a wide 
range of cosmologi cal data. Specifically , AW = 5 is within 
the allowed region (Hamannetal. 1 120 lot iDunklev et all 12010; 
iGiusarma et al.120111) . The allo wed mass range of the sterile neu- 
trinos has also been studied. In Hamann et al. (2010) a scenario 
with massless flavour neutrinos and two thermalized sterile neu- 
trinos with a common mass was investigated. Using a combina- 
tion of different cosmological probes, including the CMB and 
the galaxy power spectrum, they found a 95 % upper bound of 
0.45 eV for each sterile neutrino. For a single sterile neutrino 
the upper bound was 0.48 eV. There is some tension between 
these limits an d the best-fit masse s of the sterile neutrinos that 
we adop t from IKopp et al] (1201 ll) . However, we note that the 
limits in lHamann et alJ (1201 Oh are derived assuming a flat uni- 
verse with a cosmological constant, while we will allow for both 
spatial curvature and a dark energy equation of state, w + - 1 in 
our analysis. 

3. Method 

We study two different cosmological scenarios, the standard flat 
ACDM model, and an extended model where we allow the dark 
energy equation of state and the spatial curvature of the universe 
to vary, which we will refer to as wCDM+k. For both of these 
models we estimate the model parameters when including zero, 
one (1 v s ) or two (2 v s ) sterile neutrino species. We assume the 
sterile neutrinos to be fully thermalized, and adopt the best-fit 
m asses for the steril e neutrinos from reactor experiments found 
in iKoppet alJ (120111) . that is m v — 1.33eV in the 1 v s scenario, 



and m v ' — 0.68eV and nv v — 0.94eV in the 2 v s scenario. 
We assume the three species of flavour neutrinos to be massless, 
which should be a good approximation when the sterile neutrino 
masses are in t he high end of their co smologically allowed mass 
range (see e.g. Giu sarma et al.l (1201 ll) ). 

We use a modified version of the publicly available cos- 
mological Markov ch ain Monte Carlo sampler CosmoMC 
dLewis & Bridld 120021) to compute the parameter limits. For 
the ACDM model we vary the parameter set {a>t,, w f , 9, t, n s , 
In 10 10 As }, and for the wCDM+k model we also include w and 
Clk as free parameters. a>t> and co c are the physical baryon and 
cold dark matter densities, respectively. 9 is the ratio of the sound 
horizon to the angular diameter distance, r is the optical depth, 
n s and A s are the primordial scalar spectral index and amplitude 
(at k = 0.05Mpc~'). w denotes the dark energy equation of state 
(assumed to be constant), and is the curvature density. For 
exact parameter definitions we refer to the CosmoMC code. We 
marginalize over the SZ amplitude. All the listed parameters are 
given flat priors. 

We also use two different combinations of data sets. First, 
we only use CMB data from the WMA P 7 year data release 
( Komatsu et alJ201 lllLarson et al .1120111) . which we will refer to 
as WMAP7. Then we also include data on large scale structures 
from the Sloan Digital S ky Survey DR7 luminous red galaxy 
sample dReid et alj 120101) . Supernova 1A data from SDSS-II 
dKessler et alJ 120091) and a prior on the Hu bble parameter of 
H = 73.8+2.4km s^Mpc -1 dRiess etalJ201ll) . We willreferto 
this combination of data sets as WMAP7++. For the wCDM+fc 
model we only use the WMAP7++ data sets, as WMAP7 data 
alone have very little constraining power for this extended pa- 
rameter space. 



4. Results 

In Figure [TJ we show ID marginalized probability distributions 
for a few selected parameters for the ACDM, ACDM+lv s , and 
ACDM+2v. s cases, and the corresponding numerical limits are 
given in Table [1] 

Interestingly, we see that including sterile neutrinos with 
their masses fitted to reactor experiments will shift the age of 
the universe significantly. While ACDM favours an age of the 
universe of 13.75 + 0.13 Gyr with WMAP7, the inclusion of 2 
sterile neutrinos leads to a preferred age of only 12.77 ± 0.11 
Gyr. For the WMAP7++ data, the corresponding age estimates 
are 13.70 ± O.lOGyr and 12.55 ± 0.09 Gyr. The reason for these 
large shifts can be found in the corresponding shifts in the dark 
energy density, £2a- The shift in Qa can be understood by consid- 
ering changes in the time of matter-radiation equality, feq. At the 
time of equality the sterile neutrinos were still relativistic, thus 
they contributed to the relativistic (radiation) energy denisity, 
and keeping all other parameters constant, additional sterile neu- 
trinos will shift feq to later time s. The CMB power spectru m 
is quite sensitive to feq (see e.g. lLesgourgues & Pastorl (120061) ). 
and to shift feq back, the matter density must be increased. When 
we require a flat universe, this will lead to a reduction of £2a and 
thus a younger universe. 

In Figure [2] and Table [2] we show the corresponding results 
for the wCDM+k model, but only for the WMAP7++ data sets. 
We basically find the same shifts in Qde (corrsponding to £2a) 
and in the age of the universe, for the same reasons as explained 
above . When we apply the tight prior on Hq from Ries s et aT] 
j2"0Tlt) . WMAP data will constrain the universe to be close to 
flat for all models, as can be seen in the resulting limits on Q^. 
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Fig. 1. Marginalized parameter distributions for ACDM with and without additional sterile neutrinos. Solid, red lines: No sterile 
neutrinos. Dashed, blue lines: 1 sterile neutrino. Dotted, black lines: 2 sterile neutrinos. Thin lines show results from using WMAP7 
data only. Results from WMAP7++ are shown with thick lines. 



Parameter 


ACDM 


ACDM+lv, 


ACDM+2v s 


WMAP7 





0.970 


+0.014 


0.958 


±0.013 


0.962 


±0.013 


n A 


0.707 


+0.030 


0.646 


± 0.035 


0.598 


± 0.039 


Age (Gyr) 


13.75 


±0.13 


13.28 


±0.12 


12.77 


±0.11 


H (km/s/Mpc) 


72.0 


±2.5 


67.8 


±2.2 


68.0 


± 2.2 


A* 2 







-3.4 




-3.6 




WMAP7++ 


n s 


0.974 


±0.012 


0.982 


±0.012 


0.987 


±0.012 


Qde 


0.699 


±0.014 


0.668 


±0.015 


0.640 


±0.016 


Age (Gyr) 


13.70 


±0.10 


13.04 


±0.10 


12.55 


±0.09 


H (km/s/Mpc) 


71.5 


±1.1 


70.2 


±1.1 


71.3 


±1.1 


A* 2 







-24.7 




-22.6 





Table 1. One-dimensional marginalized parameter limits for ACDM with and without sterile neutrinos. The errors shown are the 
lcr deviations from the mean value. We also show the Ax 2 values for the different models where A^ , j = x\cdm "^model- 



When including sterile neutrinos, we notice that w is shifted 
into the w < — 1 phantom regime. For decreasing Qde, w is 
forced to smaller, i.e., more negative values to obtain the late 
time acceleration required by the supernova data. 

In Figure [3] we show 2 dimensional 68% and 95% confi- 
dence contours for a few strongly correlated parameters. We see 
from the Q^-vv contours, that for no sterile neutrinos, the ACDM 
model (Q.k = and w = -1) falls within the 68% contours, 
while with sterile neutrinos, the ACDM +v s model falls just out- 
side the 95% contours. This indicates that, if further oscillation 
experiments confirm the existence of sterile neutrinos with prop- 
erties close to what was found in iKopp et ail (1201 lh . this will 
imply some tension between the cosmological constant as dark 
energy and cosmological data. One might argue that models with 
w < -1 are unphysical. However there exist several physically 



more viable dark energy models that will give an effective equa- 
tion of state w < — 1 at late times and thus be able to accommo- 
date the current cosmological data better than the ACDM model 
in a cosmology with more mass in the neutrino sector (see e.g. 
lLa Vacca et al] (120091) ). A discussion of different dynamical dark 
energy models are outside the scope of this paper, but our re- 
su lts indicate that st erile neutrinos with the properties derived 
in Kopp e t alJ (1201 lh would call the cosmological constant into 
question as the explanation for dark energy. 

In Tables Q] and [2] we also show the Ay 2 between the best-fit 
models in the cases with and without the added massive sterile 
neutrinos. We see the addition of sterile neutrinos decreases the 
fit to the data, especially when using the WMAP7++ data sets 
and the ACDM model, giving A% 2 = -22.6 for the 2v s case. 
When opening for the £2^ and w degrees of freedom, the A% 2 re- 
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Fig. 2. Marginalized parameter distributions for wCDM+k with and without additional sterile neutrinos. Labels are the same as in 
Figure[l] Since these models are very poorly constrained when using WMAP7 data only, we only show the results from WMAP7++. 



Parameter ACDM ACDM+lv, ACDM+2v s 



n s 


0.968 


±0.014 


0.959 


±0.013 


0.963 


±0.013 


Qde 


0.700 


±0.014 


0.668 


± 0.015 


0.638 


±0.018 


Age (Gyr) 


13.87 


±0.31 


13.05 


±0.25 


12.58 


±0.26 


H Q (km/s/Mpc) 


71.7 


±1.5 


72.7 


±1.4 


73.6 


±1.9 


£l k 


-0.0033 


±0.0066 


0.0021 


±0.0062 


0.0016 


±0.0071 


w 


-1.11 


±0.10 


-1.29 


±0.12 


-1.33 


±0.14 


A* 2 







-12.0 




-9.3 





Table 2. One-dimensional marginalized parameter distributions for the wCDM+k model with and without sterile neutrinos. 
The errors shown are the lcr deviations from the mean value. We also show the Ax 2 values for the different models where 

^model _ ^ACDM ~~ ^model' 



duces to -9.3, underlining the need to look beyond a flat ACDM 
model if the existence of the sterile neutrinos in Kop p et aT] 
J20T1 is confirmed. It is interesting to notice, that the worsening 
of x 1 from introducing 2v s in our analysis, is of the same order 
as the impr ovement of the fit by introducing the same sterile 
neutrinos in Kopp et al.l d201 ll) . As previously stated, our point 
of departure in this work is that neutrino oscillation experiments 
are less prone to systematical errors than cosmolog_ical observa- 
tions, justifying the use of the results from Kop p et alj d201 lb as 
an input in the cosmological models. 

5. Discussion and conclusions 

To summarize, we have investigated how the presen ce of one or 
two st erile neutrinos with the properties estimated in Kopp et aT] 
changes the preferred values of cosmological parameters. 
Rather than deriving constraints on neutrino properties from cos- 
mology, we chose the opposite approach of using neutrino ex- 
periments to constrain cosmology. We think our approach can 
be justified, since the uncertain factors in cosmology, like the 



nature of the dark energy, are arguably larger than those in neu- 
trino physics. 

We analyzed the ACDM model and the wCDM model with 
spatial curvature as an added parameter, and for two data sets: 
WMAP7 alone, and WMAP7 plus large-scale structure, super- 
novae type la, and the HST result for the Hubble parameter. The 
most interesting changes from the standard ACDM model with 
no sterile neutrinos were in the equation of state parameter w and 
in the age of the universe. 

In the wCDM+A:+2v s model we found a preferred age of 
~ 12.5Gyr. One might question whether such a young uni- 
verse would be in conflict with other cosmological observations. 
Regarding observations of high redshift objects, which ages are 
derived from their redshifts, this should not be a problem. As 
an example, the age of the universe at redshift 12 in the mean 
parameter value wCDM+k+2v s model from WMAP7++ is only 
16% lower than for a standard ACDM model. 

One may also put lower limits on the age of the universe 
by measuring the age of the oldest objects in a cosmology- 
independent way. A common way to do this is by main sequence 
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-1.8 -1.6 -1.4 -1.2 -1 -0.8 -0.02 0.02 

w n K 

Fig. 3. Marginalized 2D parameter distributions for wCDM+k with and without additional sterile neutrinos. Labels are the same 
as in Figure [T] Since these models are very poorly constrained when using WMAP7 data only, we only show the results from 
WMAP7++. 



fitting in globular clusters. In iGrattone t al. (2003) they use this 
technique, and estimate the age of the oldest globular cluster in 
the galaxy to be 13.4±0.8±0.6Gyr (statistical/systeme tic errors), 
which leaves the 2v s models within the error bars. In Wan g et all 
(2010) nine clust er with estimated a ges of around 14 Gyr were 
found in M3 1 . In lFrebel et al.l d2007l) they use decay rates of ra- 
dioactive isotopes in nearby stars to estimate a stellar age of 13.4 
Gyr. However, the uncertainties in these kinds of measurements 
are large, and the authors estimate an uncertainty of ~ 2 Gyrs, 
which also places the model with 2 massive sterile neutrinos well 
within the allowed range. 

We found that the 3 active + 2 sterile neutrinos scenario 
prefers an equation of state parameter for dark energy w < — 1, 
with the cosmological constant being ruled out at 2<x. The fact 
that w < — 1 should not be taken as an indication of phantom 
energy. It more likely means that the correct dark energy model 
cannot be described by a constant w. If the evidence for sterile 
neutrinos from oscillation experiments becomes conclusive the 
implication could be that the cosmological constant is ruled out 
as dark energy. 

Acknowledgements. We thank Frode K. Hansen for useful discussions. The 
results described in this paper have been produced using the Titan High 
Performance Computing facilities at the University of Oslo. 
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